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Abstract 

We extend the Ito-to-Stratonovich analysis or quantum stochastic dif- 
ferential equations, introduced by Gardiner and CoUett for emission (cre- 
ation), absorption (annihilation) processes, to include scattering (conser- 
vation) processes. Working within the framework of quantum stochas- 
tic calculus, we define Stratonovich calculus as an algebraic modifica- 
tion of the Ito one and give conditions for the existence of Stratonovich 
time-ordered exponentials. We show that conversion formula for the co- 
efficients has a striking resemblance to Green's function formulae from 
standard perturbation theory. We show that the calculus conveniently 
describes the Markov limit of regular open quantum dynamical system- 
sin much the same way as in the Wong-Zakai approximation theorems of 
classical stochastic analysis. We extend previous limit results to multiple- 
dimensions with a proof that makes use of diagrammatic conventions. 



1 Introduction 

Quantum stochastic calculus was developed as a framework to construct 

concrete models of irreversible quantum dynamical systems. Prior to this, mod- 
els tended to consider a system couple to an environment, but with only the 
system being accessible to physical measurement: as a result the environment 
observables were often relegated to a secondary status, leaving one with a master 
equation for the state of the system only [Hj. 

Hudson and Parthasarathy |3 in 1984 presented a rigorous theory of inte- 
gration with respect to processes on Bosonic (later Fermionic) Fock spaces gen- 
eralizing the Ito-Doob theory of stochastic integration. In addition to integrals 
with respect to time, they also introduced integrals with respect to creation, 
annihilation and number (more generally, scattering) processes. Motivated by 
non-commutative Feynman-Kac formulae, they were able to describe unitary 
dynamical evolutions of a system coupled to the Fock space environment which 



reduced to an irreversible Markov dynamics for the system when averaged over 
the partial trace with respect to the Fock vacuum. Here the Schrodinger equa- 
tion is replaced by a quantum stochastic differential equation (QSDE) driven 
by the creation, annihilation and scattering processes. 

Independently, Gardiner and CoUett [J] in 1985 gave the version of quantum 
stochastic integration for the Bosonic creation and annihilation that is best 
known amongst the physics community. Although they did not include the 
scattering processes, they did introduce several important physical concepts, 
in particular, they gave to the noise the status of a physical observable. This 
has been vital in subsequent analysis of quantum filtering and feedback where 
the environment can act as an apparatus/ communication channel |7]-|ni- In 
their analysis, they also introduced the Stratonovich version of the theory by 
extending the usual mid-point definition to non-commuting processes. This 
is a natural physical choice for two reasons: unlike the Ito form, the Leibniz 
rule of differential calculus holds for Stratonovich differentials and so physical 
symmetries are more apparent; secondly, in classical analysis it is generally the 
case that if the model can be obtained as a singular limit of regular dynamical 
models, then it is the Stratonovich form that resembles the pre-limit equations 
the most. Historically, it was actually the latter reason that lead to Stratonovich 
initially introducing his modification of the Ito theory. In ordinary stochastic 
analysis, such results involving a central limit effect for stochastic processes are 
known as Wong-Zakai theorems. Our original motivation stems from quantum 
Markov limits - ^S] and the desire to understand the limit processes in a 
Stratonovich sense. 

To give a concrete mathematical account, we start from the quantum Ito the- 
ory developed by Hudson and Parthasarathy, and deduce quantum Stratonovich 
calculus as an algebraic modification of the quantum Ito one which restores the 
Leibniz rule. An approach starting from Gardiner and Collett's input processes 
would have been more appealing from a physical point of view, however, we 
do not want to bypass questions of the mathematical status of the objects con- 
sidered. Traditionally, the Stratonovich integral is defined through a midpoint 
Riemann sum approximation: this has been extended by Chebotarev |12 | to 
quantum stochastic integrals, however, we emphasize that our formulation here 
is to define Stratonovich integrals as combinations of well-defined Ito integrals 
where possible. Rather than there been an unique version, we find that there 
are degrees of freedom in how we actually achieve this - we refer to this as 
a "gauge" freedom - and that the standard (symmetric) choice, corresponding 
to the midpoint rule, is just one possibility. We give the self-consistency for- 
mula G = Go -I- GqVG relating the matrix of Ito coefficients G to the matrix of 
Stratonovich coefficients Gq. (Here the "potential" V is half the noise covariance 
matrix plus the gauge.) Rather surprisingly, this has the same algebraic form as 
the one relating the free and perturbed Green's functions in scattering theory: 
a fact that we readily exploit. It is shown that the Ito coefficients of a unitary 
process are related to Hamiltonian Stratonovich coefficients ^GJ — — Gq^ and 
that this is true for any gauge so long as the self-consistency formula can be 
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solved. This allows the interpretation that the Stratonovich calculus can be 
viewed as a perturbation of the Ito calculus, and vice versa. 

A major motivation here is the results of |15j for a quantum Markov limit in- 
volving emission, absorption and scattering. We formulate the limit as a Wong- 
Zakai result where the Stratonovich QSDE resembles the pre-limit Schrodinger 
equation (with gauge set by the imaginary part of the complex damping) . Our 
key requirement is convergence of the "Neumann series" G = X)n=o (^CIq)". 
We generalize the result to multiple channel noise sources and make the proof 
more accessible by means of diagrammatic conventions which make the connec- 
tions with the Dyson series expansion transparent. We also present the results 
in a fluxion notation, as an alternative to the differential increment language, 
which is closer to the formulation employed by Gardiner and CoUett Q and 
effectively generalizes their results to include scattering. This also reveals a new 
representation (lemma 3) for the Evans-Hudson flow maps. 



2 Quantum Stochastic Calculus 

2.1 Quantum Processes 

Given a Hilbert space the Fock space over f)i is the Hilbert space F (l)i) 
spanned by symmetrized n-particle vectors Lpi^ ■ ■ ■ '^(p„ := (n!) X^o-es '/'o-(i)® 
• • • (g) <Po.(n) for n > arbitrary, (p^ G f)i and 6„ the group of permutations 
on n labels. The special case n = requires the introduction of unit vector 
Q called the Fock vacuum vector. The inner product on F([)i) is given by 

The exponential vector with test function (p E i)i is defined to the Fock space 
vector 

oo ^ 

= XI — y ^ ■ ■ ■ ^ (1) 

n told 

and we have {s {(p) \s{ip)) = exp{ip\ilj). If 5 is a dense subset of f)i then the 
vectors e {ip), with ip E S, are total in F (f}i), that is, the closure of the span of 
these vectors gives the whole Fock space. 

We now take the one-particle space to be the Hilbert space of C^-valued 
square-integrable functions of positive time: this consists of measurable func- 
tions f = (/i,-- - ,/„) with /o°° I]r=i l/i (^) P < oo- All orthogonal projec- 
tion Hs is defined on the one-particle space for each s > by taking H^f = 

(x[os]/ir'' 7X[os]/rt) where X[a s] the indicator function for the interval 
[0, s]. An A^-channel quantum noise source is modelled by operators processes 
I : i > o| acting on the corresponding Fock space ^. For a,(3 £ {0, 1, • • • , n}, 
these processes are defined on the domain of exponential vectors by 



e(f) 



(^Af-l^f:{s)gp{s)ds'^sig)'^=0 (2) 
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where we include the index zero by setting /o = ffo = 1- We shall adopt the 
convention that lower case Latin indices (with the exception of t and s which we 
reserve for time!) range over the values 1, • ■ • , iV while lower case Greek indices 
range over 0, 1, • • • , iV. We also apply an Einstein summation convention for 
repeated indices over the appropriate range. 

We also fix a Hilbert space [}, called the initial space. Let D he a, domain 
in i) and take S to be the space of bounded C-'^-valued functions. A family 
X. = {Xt : t>0} of operators on () (g) 5" is said to be an adapted quantum 
stochastic process based on {D, S) if, for each t > 0, Xtu (g) e (f ) is defined for 
each u G D and f = (/i, • • • , /„) £ S and is independent of the values fi (s) for 
s > t. If X"^ are adapted processes, their stochastic integral X. may be written 
as (implied summation!) Xt = aj, (s) X"f^ap (s) ds, with the meaning that 



M(g)e(f) 



ai (s) Xfap (s) dsv®e (g) 

/: (s) {u ® e (f) \xfv ® £ (g) ) gp [s) ds 



The stochastic integral is more often written in the form Xt — X"^dA"^ 

or equivalently dA"^ X"^ with these integrals making sense as Riemann- 
Ito limits for locally square-integrable integrands. The stochastic integral will 
again be an adapted process. At the moment, the symbols aj, (t) , a^g (t) have no 
meaning other than notational, however aj (t) and oq (t) are easily interpreted 
as the identity operator and Ui {t) as a pointwise Malliavin gradient. What is 
crucial is that they appear in Wick order: that is, aj, (t) to the left of a/3 {t). Let 
Yt — Jq ojj (s) Y^^ap (s) ds be a second integral, with the Y""^ adapted, then 
we have the formula [3] 

XtYt = / ai{s)X'^^Ysap{s)ds+ [ ai{s) X^Y^f" ap {s) ds 
Jo Jo 

{s)XrP,,Yfap {s)ds, 







where we introduce 



P^^- = l ^' f'-"^^' (3) 
' 0, otherwise. ^ ' 



Introducing the differential notation dXt — X"^dA"^ , etc., we may write 
the quantum Ito formula in the more familiar guise as 

d {XtYt) = {dXt) Yt + Xt (dYt) + (dXt) (dYt) (4) 

where the Ito correction is (dXt) (dYt) = X^'Yt'^dAt'^ = Xt'" P'^'^Yt'^ dA"'^ . 
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Let us define iterated integrals in the natural way: 
(£(f)| / dAtf----dAl^^^ £(g)) 

= f fL (in) •••/:, ih) 5/3„ (tn) • • • ffft ih) {e (f) I 6 (g)) (5) 

where A„ {t) is the simplex t > i„ > • • • > ii > 0. 
2.2 Quantum Markov Evolutions 

Let {Gap} be bounded operators on a fixed initial space [), then there ex- 
ist a unique solution U. to the equation [/t = 1 + J* GapUsdA^f^ which we 
can naturally interpret as the QSDE dUt = GapUtdA"^ , with Uq — 1. In 
such cases, we may write U. as the Dyson-Ito time-ordered exponential Ut = 
fioe^p S^GapdA^^. 

Proposition 1: Let U. be the solution to the QSDE dUt — GapUtdA^^ , with 
Uq — 1, where the {Gafi} are bounded operators on (). Necessary and sufficient 
conditions for U. to be unitary process are that Gap + G^^ + G\^Gip = = 
Gap + Cjjj^ + GiaG\p with the general solution 

Gij = — 5ij; Gio — Li; 

Goj = -LiWkf, Goo^~^LlLk-iH. 

where Wij , Li and H are bounded operators on the initial space with H self- 
adjoint and W^jWjk = Sik = WijWj^. 

Proposition 2: Let U. be the unitary process described above The corresponding 
flow map is given by Jt (X) :— Uj {X (g) 1) Ut for bounded X on the initial space. 
We find that Jt (X) satisfies the QSDE 

dJt{X)^Jt{Cap {X))dAf, 

where the Evans-Hudson maps '16^ are given by Cap {X) :— XGap + G'^^^X + 

Gla^Gip. 

Following Lindblad ^7] , the dissipation of a linear map C on the algebra of 
bounded operators on f) is defined to be the bilinear mapping DC : (X, Y) n- 
C [XY] - C{X)Y - XC [Y). 

Proposition 3: The Evans-Hudson maps for a unitary flow satisfy Cap {X)^ = 
Cpa (^^) o.nd their dissipation is described by the equation 

DCaP {X, Y) = Ca^ (X) P'^^'Cp (Y) . 
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2.3 Approximations 

For each A > 0, wc set {t, A) = and take af (t, A) = {ip (i, t, A)) where 
(fi {i, t, A) is a C-'^-valued square- integrable function on [0, oo): we also take t 
(f (i, t, A) to be strongly differentiable. We assume that (ip {i, t, A) \(p (j, s, A)) = 
Cij {t — s, A) where Cij (t, A) = Cji (— t) is integrable in r, and that we have the 
convergence liniA^o (''': ^) = ^ij ^ {t) in the sense of Schwartz distributions. 
We set Kij = liniA^o /q°° Cij {t, A) dr and note the identity Kij + Kj^ = Stj . 

The af {t, A) are approximations to quantum white noises. We may intro- 
duce integrated processes A"^ (A) := aj^, (s, A) (s, A) ds which serve as ap- 
proximations to the fundamental processes. The approximation will be termed 
symmetric if Kij = ^Sij, but this is only a special case. Defining smeared 
exponential vectors by 

£A(f)=exp|^ J2fi{t)al{t)dtjn, (6) 

we see that the limit 

lini (£a (f) I / a* {tn, A) • • • a* (ti, A) 0/3 A) • • • ap^ {ti,X) ex (g)) 
-^'^o JA„(t) 

coincides with (5). Therefore, whenever we consider the limit of Wick ordered 
expressions, it doesn't matter whether we have the symmetric approximation or 
not. 

A term such as /^^(t) ^» (^2, A) a] (ti. A) must, however, be put to Wick order 
as /^^(t) at {ti,X) tti {t2, A) plus an additional term /^^(^) Cjj (t2 — *2) A) which 
converges to k^ L As a rule, expressions out of Wick order will have a limit that 
depends on the constants K^j. 



2.4 Notation and Conventions 

Let X = {Xap} be an (A^ -|- 1) x (A^ + 1) matrix of bounded operators on some 
Hilbcrt space. (Previously, we had the Ito coefficients which were operators on 
f) tJ.) We shall adopt the matrix representation 



X = 



/ A, 



00 



A'lo 
A'20 



Aoi,Ao2,--- \ 



X 



(7) 



where X is the N x N sub-matrix consisting of the entries {Xij}. 
For instance, P = {P"''} will be a projection operator: 















, Q 





^ 





^ 



(8) 
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The quantum Ito correction is therefore described by the matrix X (t) PY (t) = 
{X^^ (t) P^^'^Y^p {t)}. 

Let {Gafi\ = G be a matrix with bounded operators on [) as entries, then 

the Dyson-Ito time-ordered exponential Ut = T^di |exp GapdA°'f^^ will be 
unitarity if, from proposition 1, 

G + Gt + G^PG = 0, G + Gt + GPG^ = 0. (9) 

It is relatively easy to see that the Ito coefficients then take the general form 



PGP 
QGP 



W P, 
-(PGQ)tw, 



QGQ = i (PGQ)^ (PGQ) i 



H 












(10) 



where W'I'W = P = WW^ (i.e., the restriction of W to () (g) is unitary) and 
H is self-adjoint on (). 



More explicitly, we may set 



W 







PGQ 



L2 



V 



o,o,--- \ 



where Wij, Li and H are the operators on t) introduced in proposition 1. 

We should remark that the restriction to a finite number N of channels is 
not essential and that the unitary process exists under certain conditions on the 
boundedness of G as a matrix operator [2] . 



3 Quantum Stratonovich Calculus 

We wish to write the quantum Ito formula in the form 

d {XtYt) = {dXt) oYt+Xto (dYt) . (11) 
This can be achieved by formally defining 

{dXt)oYt : ^idXt)Yt+X^^VrYt''^dAf, 

Xto{dYt) : =XtidYt)+Xnvr)^Y,'"'dAf, (12) 

where the |Vj"''| may in general be taken as adapted processes, however, we 

shall take them to be just scalar coefficients. It follows that [{dXt) o Yt]^ = 
Y^^odX}, and we recover the Ito formula provided we have the condition -f 
(yu^y ^ p^iu rpj^g simplest solution possible is to take V^"" = ^Pf"" and 
this corresponds algebraically to the traditional Stratonovich definition of a 
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differential. The general solution however takes the form V^" = ^P^" + iZ^", 
where the constants {Z^^"} satisfy [Z^^")* = Z"^. The appearance of these 
constants is similar to the ambiguity in the Tomita-Takesaki theory, and we 
refer to them as a gauge freedom. We shall identify the F*-' with the constants 
Kij occurring in the approximation scheme. Let us take V = {F"''} to be the 
family of constants, then the requirement is V + = P with general solution 

V = + iZ where = Z. We shall take the [Z"^] to be scalar constants 
and set Z^° = Z''° = Z°^ = 0. This implies that Z = ZP = PZ and so 

V = VP = PV : 












V 



where V = ^ + iZ with = Z. Note that V is a normal operator with VV^ = 

VtV = 4-i + Z2. 

It should be pointed out that we have the relation 

dJt {XY) = {dJt (X)) o Jt {Y) + Jt (X) o {dJt (F)) 

which we can get either from taking differentials of the homomorphic prop- 
erty Jt (XY) = Jt {X) Jt (Y), or explicitly by noting that {dJt {X)) o Jt (Y) = 
J((£a/3 {X) Y + £ai^ {X) V^Cnij {Y))dA^^ , etc, and using proposition 3. 



3.1 Stratonovich-Dyson Time Ordered Exponentials 

Let us now suppose that U. is simultaneously the solution to the Ito QSDE 
dU = {dG) U, with dG = GafjdA""^ as before, and a Stratonovich QSDE (for a 
fixed gauge!) 

dUt = {dGo{t))oUt, Uo = l, (13) 
with dGo = GlfjdA'^f^. In such cases, wc may write U. as the Dyson- Stratonovich 
time-ordered exponential Ut — Tsd |cxp/o ^a£('^^"^| ^^"^ shall refer to Gq = 
l^a/s} ^^'^ matrix of Stratonovich coefficients. 

Self-consistency requires that dU — (dG'^) o U = (dG) U and so we should 
have that dU = (^6°) U+G'^^V'^G,^ UdA<^^ = (g^^ + G^^V^^G^p) U dA<^0 . 
This means that the Ito coefficients G = {Gajs} are related to the Stratonovich 
coefficients Go = |g°^| by 

G = Go GoVG. (14) 

As we shall see, so long as 1 -|- VGP is invertible, we may solve for Go in 
terms of G. Similarly, invertibility of 1 — PGoV implies that we may write G 
in terms of Go. It might be remarked that the relation (14) also applies if we 
consider matrices G, Go of adapted processes. 
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What is rather astonishing is that relation (14) is precisely of the form 
relating free and perturbed Green's functions. Let us recall briefly that if 
H = Hq + F is a Hamiltonian considered as a perturbation of the free Hamil- 
tonian Hq then the resolvent operator Q [z) = (z — H)~^ is related to the free 
resolvent (z) = {z — Hq)~^ by the algebraic identity 

= 00 + 0oV0 (15) 

for all z outside of the spectra of H and Hq. The identity may be rewritten as 
= {I — 0oy)~^ 00 and iterated to give the formal expansion — 0o + 0oy0o + 
0oy0oy0o + ■ ■ ■ which, when convergent, is the Neumann series. The details 
of the actually scattering are contained in the operator T := V + V0V and we 
have the identity = 0o + 0o'^0o- 

3.2 The T-matrix 

We now exploit the similarity between (14) and (15). We begin by introducing 
the operator 



V + VGV 















(16) 



Assuming that 1 — VGqP is again invertible, we obtain the following identities 

GoT = GV, (18) 

TGo = VG, (19) 

G = Go + GVGo, (20) 

G = Go + GoTGo. (21) 

The proof of (17) comes from writing 

T = V + T (Go + GoTG) V = V + VGqT 

so that (1 — VGqP) T = V. The remaining identities are just precise analogues 
of well-known relations for resolvent operators . 

Combining (17) and (21) we see that G can be expressed in terms of Go as 

G = Go + Go(l-PVGoP)-'VGo. (22) 

In particular, we see that G is bounded. We may then invert to get 

Go=G G(1 + PVGP)^VG. (23) 

The equations (22) and (23) reveal a remarkable duality between the Ito and 
Stratonovich coefficients. (Of course this just means that we may view either 
as a "perturbation" of the other!) 
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If PVGoP is a strict contraction, then we may develop a Neumann series 
expansion G = Go + GqVGo + GqVGoVGo + • • • = Gq EZo (VGq)". 
It is convenient to introduce a related matrix 



TGo 



VG 



(24) 



so that G = GoF. 

3.3 An "Optical Theorem" 

Let us next suppose that the Stratonovich coefhcients take the Hamiltonian 
form 

Go= - iE 



where E is a bounded, self-adjoint operator on f) 



(25) 

We then have the 



relation G, 



-Go and set 



PEP 















so that E is self-adjoint on [) (g) C^. When our invertibility condition is met, it 
is easy to see that matrix T exists and can be written as 



T=[l + iVE]-W^ — 



1 



iE 



(26) 



(In the special case where Z — 0, the self-adjointness of E ensures that 
V = l-l-i^E is invertible by von Neumann's theorem '53'. Therefore the existence 
of matrix T is guaranteed. More generally, so long as the value 1 lies in the 
resolvent set of ZE, this theorem implies the existence of T.) 

The related matrix F then takes the form 

/. 



1 


0,0, ••• 






-iT2jEjo 


LL 



with F = 1 - iTE = TV-i, so that F = [1 -f^ iVE]"\ 

Lemma 1 ("Optical Theorem"): RcT > and in particular T satisfies the 
identity 

T-hT^ = FFt = F^F. (27) 
Proof. We have thatT -I- may be written as 



FV-l-VtFt = F 



V (l - iEVt) + (1 + iVE) V^j Ft = F [V + V^] F^ = FF^ 
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It is then relatively straightforward to show that 

- (l-iVtEHl + iEV) -^^-^^-^^ + ^^(^^)^] 

^FtF. 



(l + iEV)(l-iVtE) 
■ 

A similar calculation shows that ImT = — TET^. 
3.4 Unitarity 

We now wish to show that the choice of Hamiltonian Stratonovich coefficients 
naturally leads to unitary processes. 

Lemma 2: Let Gj = — Gq be bounded with 1 — VGqP invertible, and set 
G° (t) = J*GlpdA°''^. Then the solution to the Stratonovich QSDE dU = 
(dG*^) o U , Uq — 1 will be unitary. 

Proof. This fact is an immediate consequence of the optical theorem (27). 
To establish the isometric property for the Ito coefficients, first observe that 
G + Gt = -Go (T + Tt) Go while 



GPGt = GoFVFtGj, = Gq 












FF^ 



gJ, --Go (T + Tt)Go. 



The isometry condition in (9) then follows from the first part of (27). The 
co-isometric property likewise follows from the second part. ■ 



3.5 Changing Gauge 

Let G be a fixed Ito coefficient matrix related to the Stratonovich coefficient ma- 
trices Gq"'' and Gq"'' with gauges Z^"^ and Z^^^ respectively. The two matrices 
will then be related by the perturbative formula 

g'^;^ = gI,'^^ + iG[,"^ (zt-^) - Z'') g[,°\ 

In particular, we can relate Gq"' for a non-zero gauge to the symmetric (gauge 
zero) form Gq"'. As we shall see, the gauge Z has the physically interpretation 
as the imaginary part in the complex damping and in many applications this 
may be small |2(J| . 



4 Wick Ordering Rule 

Let X. and Y. be quantum stochastic integrals with adapted integrands as be- 
fore. In terms of our notation involving the aj^, (t) , ap (t), the product XtYt = 
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/g dsi /q ds2 ajj (si) X^fa^ (si) aj^ (52) ^t-^cLv (S2) is not immediately interpreted 
as an iterated integral since it is out of Wick order. However, the rule for achiev- 
ing this is formally equivalent to the kinematic relations 

po^t^Yl'^a^, {t), t< s; 
[a„ (t) ,Ys]^{ V^^^yra, (t), t = s; (28) 
0, t>s: 



under the integral sign, along with its adjoint 



Xt,al{t)'] - [ai3it),Xl] \ These 
da (t) , aj; (s) 



relations can be viewed as the formal commutation relations 

¥"^6+ {t - s) + {V^'^y S_ [t - s) at work, where the 5± are one-sided delta- 
functions: / 5± (/) f (t) = f (0^), see e.g. It is possible to interpret the 
(0 quantum white noise operators, but we do not stress this point further 
here. At this stage, we could switch to a fluxion notation such as Xt — = 
a]j (t) X^^ap {t), etc., and write the Ito formula as ^ {XtYt) — Xt o Yt + Xt o Yt 
with the convention that 



XtoYt^ ai it) Xfap {t) Yt = 4 {t) xfYtap (t) + at {t) X^V^^Yi^^ap {t) . 

In particular, we have the following interpretation of the results of the pre- 
vious section: The equation Ut = Go {t) Ut with Gq {t) — aj^ (t) G^^pap {t) is out 
of Wick order, but can be put to Wick order as Ut = aj, [t] GapUtap (t). The 
relation [a„ (t) , Ut] = V'^G'^'Uta^ (t) imphes that 

a„ (t) Ut = {da^p + V^G^.p) Utap (t) = F^pUtap (t) , (29) 

where {Fap} are the components of the matrix F introduced in (24). 

The QSDE for the unitary [/., with G^p = ~iEap, will then be 

Ut = -ial (t) Eo^pap (t) Ut = -i4 (t) E^^pFp^Uta, (t) (30) 
and likewise the QSDE for the flow will be 

j^JtiX) - uHx)Ut + u}{x)Ut 

= -iU^al (t) [X, Ea^p] ap (t) Ut 

= -\al{t)ulFi^[X,E^p]Fp,Uta,{t). (31) 

Comparison with proposition 2 suggest that £^1, {X) — [X,Eap\Fp^. As 
this is an entirely new relation, we give an independent derivation in appendix 
A using only the Ito calculus. 

Lemma 3: Under the conventions and notations of the previous sections, the 
Evans-Hudson maps take the form 

C^p{X) = -iFl^[X,E^,]F,p. (32) 
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5 Quantum Wong-Zakai Theorem 



The following is the multi-dimensional version of a result first established in 

Theorem: Let a^{t,X), a = 0, 1, • • • , TV, be continuous in t creation / an- 
nihilation fields for each A > with A^^ (A) — aj^ (s, A) (s, A) ds approx- 
imating fundamental quantum stochastic processes with internal space as 
before, with fixed gauge matrix V = {F"''}. // T^'*''' = Eaf} ® aj^ (t, A) ap [t, A) 
with e'^p = Epa bounded operators on a fixed Hilbert space \) such that VE 
is a strict contraction, then the unitary family U^^^ and the Heisenberg dy- 
namical map J^^^ (X) — U^^"'^ Xu!'^^ determined by the Schrddinger equation 
tJi^^ — —iT'"^^^u[^\ t^o^'' — 1; converge in the sense of weak matrix lim- 
its to the unitary quantum stochastic process U. and corresponding quantum 
stochastic flow J. {X). The limit process U. is unitary adapted and satisfies the 
Stratonovich QSDE dUt = —\EapUt o dA"^ , Uq = 1, with gauge determined by 
V = {F"^}. 

The condition ||VE|| < 1 gives convergence of the Neumann series. It also 
implies that 1 + VEP will be invertible and therefore the Stratonovich QSDE 
makes sense. We will sketch the proof of this theorem in the Appendix B. 
Provided that the strict contractivity conditions hold, we could replace E, and 
indeed V, by suitably continuous adapted processes. 

It might be remarked that there exists an analogue of this result using Fermi 
fields in place of Bose fields • The limit QSDE changes insofar as the noises 
must now be Fermionic processes, however, the coefficients are exactly as before. 

5.1 Examples 

5.1.1 Classical Wong Zakai Theorem 

As a very special example of theorem, let us take the 1-dimensional case with 
the pre-limit Hamiltonian T^^^ determined by Eqq = H, Eqi — Eiq = R and 
Ell — with K — ^. Then the limit flow is characterized by the maps Cn {X) = 
0, Cio {X) = Cai {X) = -i [X, R] and £00 {X) = -i [X, H] - ^ [[X, R],R]. For 
the choices H = ^ {pv {q) + v (q) p) and R = ^ {pa {q) + a (q) p) , where v {■) 
and a (•) are Lipschitz continuous with \v {x) |, \a {x) \ < C (1 + |x|) for some 
constant < C < 00, and q and p are canonical position and momentum 
observables, we have that qt Jt (q) satisfies the essentially classical SDE dqt = 
[v (qt) + o- (qt) a' (qt)] dt-\-a{qt) dQt. or equivalently dqt = v (qt) dt + a {qt)odQt 
where Qt = A}'^ + A^^ is a copy of the Wiener process. 

As a result, the theorem reduces to a classical Wong-Zakai approximation 
theorem which states that, since q['^'^ = (a^ (s, A) -I- a (s, A)) ds is an essen- 
tially classical stochastic process that is differentiable in time t and converges 
almost always uniformly on compact time- intervals to a Wiener process Qt, the 
solution to the random ODE ^ = v (^X^t^\t^ + a (^X^t^\t^ q[^\ X^^^ = xa 
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similarly converges to the diffusion process X. satisfying the Stratonovich SDE 
dXt = V {Xu t)dt + a {Xt,t) o dQt, Xo = xq. 

5.1.2 Quantum Diffusions 

Taking T^-^' = a'^ {t,\) + R'^ ® a (t, X) + H leads to the limit QSDE 

dUt = -i (i? ® dAl +R^ dAt oUt 

= -i(^R(g) dAl + dAt +H^Ut- kR'^ RUtdt. 

Note that Re K = i so that dJt {X) = -Ut {[X, R]) dAl - iJt ( [X, i?^] ) dAt + 

Jt{C{X))dt where we set At = A?!, AI = A\^ and C{X) = Coa{X) = 
-i[X,iJ'] + i [R^,X] R+^R^ [X,R]. The new operator H' is H + lin{K} R^ R 
which includes an energy shift coming from the complex damping k. The theo- 
rem then reduces to a long line of results dealing with broadband noise limits, 
weak coupling limits, etc., in quantum physics. 

It is straightforward to extend this to describe coherent states, thermal states 
and squeezed states 



5.1.3 Counting Processes 

Let us consider the choice T^^'' = E ig) {a {t, X) + f (t)}^ {a (t, X) + f (t)}. This 
is what we could consider in the vacuum as an equivalent of studying E (g) 
a {t, A)^ a {t, X) in a coherent of intensity /. Here we have E^p — Ef* (t) fp {t) 
and the one dimensional form of the theorem yields the Ito coefficients Gap — 
-iEap-KEaij^^^Eip, which in this case reduce to G^p {t) = i^f^/a (t) fp {t). 
The Ito form of the limit QSDE is then 

1 + ikE 

and we introduce Nt (/) = /g (dA," + / (s) dAl" + f* («) dAf + \f{s) pds) which 
is essentially classical and corresponds to a time-inhomogeneous Poisson process 
with instantaneous rate of change ly (t) — \ f (t) p. 

Appendix A (Proof of Lemma 3) 

We restrict to the = 1 case, as this is already somewhat involved. The 
multi-dimensional case does not present any more technical difficulties. Here 
the Evans-Hudson maps are given by 

£11 (X) = W^'XW-X] 

Clo {X) = [X, L] ; £01 {X) = - [X, W; 

CooiX) = ^[L\X]L + U^X,L]-i[X,H]. (33) 
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where the operators W, H, L have the expUcit forms 

1 



w 1 - T 



--E'lO) H = Eqq + Eqi Im 



1 + ikEii 



The components of the matrix F can also be written out in detail: 
1 „ . 1 



1 + ikE II 



Fw 



1 + ikE II 



Eio, Foi = 0, Foo = 1. 



Eio- 
(34) 

(35) 



We now check that the relation (32) is correct by direct substitution. 
For a = /? = 1, we find after a little algebra that 



Cii {X) 



l + inEn ^^ l-iK*En ^ 



1 - iK*Eii 1 + ikEii 



1 - iK*Eii ' ' 1 + ikEii 
= -iFl,[X,E^,]F,,. 

For a = 1, /3 = 0, we have 



jOw (X) = 



1 + inEn 
1 - iK*En 



X,-i- 



1 



1 + inEn 



Eio 



and to compute this we need the observation that 



1 



1 + ikEu 



1 



= —IK 



1 + 'ikEii 



[X,E, 



1 



1 + IkEii 



to write 



£io {X) = IK 



[X,Ei 



1 - iK*Eii ' ' ^ 1 + ikEu 1 - iK*En 



[X,E 



lOj 



As we have £ai3 {X)^ = C^a (-'^^)) this gives the £oi {X) result also. 
The final case is the Lindbladian map Cqq. Substituting in gives 



£oo {X) = -i [X, Eoo] - i [X, Eoi Im {T} E,o] 



E, 



01 



1 - iK*Eii 



.X 



1 



1 + \kEi_i 



-E 



10 



2 1 - \K*E\i 



X, 



1 



(36) 



1 + ikEh 



-E- 
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and again we use (36) along with the observation that 



ReT= i 



1 



1 



2 1 - m*Eii 1 + \kEii 



to obtain 



£00 {X) = -i [X, Eoo] - [X, £;oi] Ti;io + EoiT^ [X, Eio] 
-lEm [X, ImT] Eio-iEoiT^ [X, En] Re {T} Ew-iEqi Re {T^} [X, En] TEiq. 
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The first three terms are -iF^p [X, Eap] F/^o, for (a, /3) = (0, 0) , (0, 1) and (1, 0) 
respectively. To tidy up the last term, we note the identities 

[X, T] ^ -iT [X, Ell] T, [X, rt] = +irt [X, En] 

so that 

[X,Imr] [X,Eii]ReT + ReT [X,En]T 
= -ir [X, En] T - [X, En] T + T^ [X, En] Re T + Re T [X, En] T 
= T^X,En]T 

and therefore the last term is -iEqiT'' [X, En] TEw = -ii^io ^w] Fw- This 
gives the desired form. 

Appendix B (Proof of the Theorem) 

We now sketch briefly the proof of limit using diagrams. Essentially, this 
is a form of the van Kampen cumulant expansion |2()j which can be described 
explicitly. The Heisenberg evolution limit is similar though more involved |15j . 

Step 1: Wick Ordering the Dyson Series 

For finite A, the mapping t i— > Ui^^ is differentiable, and we may formally 
expand as the Dyson series U^^'^ — X)n>o /a (t) ^t^'' ' ' ' ^tf^- Here A„ {t) 
is the simplex consisting of multi-times (i„, • • • ,ti) with t>tn>--->ti>Q. 

It is convenient to put the Dyson series to Wick order using the commutation 
relations. The most convenient way to describe this is to expand in terms of 
diagrammatic series and we borrow some standard techniques from field theory. 
To this end, we introduce four vertices corresponding to the four components 

of t|^^ 

2^ C. 

Scattering Emission Absorption Neutral 

The Wick-ordered Dyson series is then given by the sum 

C/f)=5]l),(A) (37) 

which we now describe. We sum over all diagrams J) obtained by writing n 
vertices in a line (?i = 0, 1, 2, • • • ) as below 



tn tj t2 ti 
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taking each of the vertices to be one of the four shown above, connecting up 
some of the creation/annihilation pairs and leaving the rest as external lines. 
With each such diagram J) we associate the operator 

2)t (A) = (-i)" £;„ • • • Si ® / Hflt (A) nC (A) Ha (A) (38) 
where, depending on the type of vertex at tj , 




Eij, scattering, 
EiQ, emission, 

Eqj, absorption, 
Eqq, neutral. 



and in the simplex integral Ila^ (A) is the product over all external lines going 
out, Ha (A) is the product over all external lines coming in and IIC (A) is a 
product over all contraction pairs. 

For instance the nine vertex diagram below 




tg ts tj Iq is ts t2 ti 



corresponds to the operator 

—iEojg Eojg Ei^oEi^jg Eoj^ Ei^oEoj^ Ei^oEoo 
«> / 4^ (t2. A) aj, (^5, A) (tg. A) 

^jgi? ~ ^7, A) Cjgig (^8 — te, A) Cjgi^ {te — t4, A) . 

The diagrams we are considering are Goldstone diagrams, or time-ordered 
Feynman diagrams j28| . The vertices however carry an operator weight E^f} 
and, even in the absence of external lines, the diagrams will not generally be 
scalars! 

Step 2: The Markov Limit of Individual Diagrams 

Let us consider the effect of the Markov limit A ^ on individual diagrams. 
If we have a contraction 

£1:^ 

over consecutive times ij+i > tj, then we will have tj+i — tj — s- 0+ in the limit. 
The effect of each such contraction is to reduce the order of the simplex by one 
degree and to introduce a multiplication factor V^^ . (Remember, we only get 
a partial contribution from the delta-function since tj+i > tj.) On the other 
hand, if any of the contracted time pairs arc not time consecutive then the 



17 



contribution converges weakly to zero! The reason is essentially that not just 
the emission and absorption times, but all the intermediate times are forced to 
be equal, and we get a collapse of the integral to zero. 

We therefore say that a diagram is time-consecutive (TC) if each contraction 
appearing is between time consecutive vertices only. 

Step 3: The Vacuum Limit 

We now &x u,v e i) and investigate the limit ]imx^o (u iSi i^\Uf^^ \v iSi ^^'^ 
which we denote by the single vertex 




We now argue that this limit will consist of the TC diagrams having no external 
lines which we can express as 

O = + • +(■■■■• •■■■■+--^--) 

+( + C:\.... + ..../Ti + ..../ZYZi.... )+ • • • 



The diagrams have been grouped by vertex number, however, it also possible 

to group them by effective vertex number which is the actual number minus the 
number of contractions and which gives reduced degree of the simplex. The 
series can be partially re-summed as 

o = + □ + + + - 

where each box is a sum over all effective one- vertex contributions: 



which is analogous to the expression of the self-energy in quantum field the- 
ory: as a sum over irreducible terms. (Note that series terminates at second 
order when there is no scattering: as this is a form of cumulant expansion, the 
emission/absorption problem is then Gaussian while allowing scattering means 
that we have cumulant moments to all orders!) Explicitly, the box at time vertex 
t corresponds to the sum 

-iEoodt + i-if EonEi.odt + {-if EquV'^'' Ei.i.V'^'^ Ei.odt + ■■■ 

= (^-iEoo - EoiV'' (rnvE) ^^°) = '^^ 

and we have summed the geometric series, convergent since |1VE|| < 1, to get 
the required coefficient Goo- We can write the expansion in the recursive form 
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+ 



DO- 



and this is interpreted as the integro- differential equation 

/ (u(g)12|dC/(Ji)(g)f^) = (w(g)i^|w(g)f^) -i / dt2lu®VL\GQQ I dUti\v®VL 
Jo Jo \ Jo 

with decaying exponential solution 

{u^n\U {t) \v^n) = (ule-""^™ It;) . 

The interpretation is intended to suggest that there is a limit object Ut 
which wc interpret as a unitary quantum stochastic process on a noise space 
with initial space f). 

Step 4: Limit of Exponential Vector Matrix Elements 

Calculating limA^o £a (f) \ U {t, A) |t; (g) (g)) docs not require too much 
machinery beyond that used in the vacuum case. Indeed, we can convert the 
new matrix elements to vacuum ones using the relation 



u®ex {t)\ui^^\v®ex (g)) = {u ® niU^-" \v ® ft 



where uj;^^ is the solution to the Schrodinger equation determined by the 
"Hamiltonian" t[^^ obtained by replacing the fields a {t, A) and {t, A) with 



at {t, A) = ai {t, ^) + J C'y (t - s, A) gj (s) , 
a]{t,X) = 4{t,X)+ j Cij{t-s)fj{t) . 

The new interaction Tp-* isn't necessarily sclf-adjoint, however, that doesn't 
effect things. We may rearrange Tj^'' in terms of the original fields as Eajs ® 
{t, A) ap {t, A) where, suppressing the A and t dependences, the are the 



operators 






Eij = 


E- 




Eio = 


Eio + Eij 


Eoj = 


Eoj + Ejk 


Eoo = 


Eqo + 


Eio 




+Ejk 


[/' 



j Cki {t-s)fi{t) 

Cu {t - s) fi (t) 



+ Eoj / Cji {t - s, A) gi (s) 



J Cu {t - s) fi (i) J Cjm {t - s, A) gm is) . 
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Therefore we only have to repeat our previous argument, but with the orig- 
inal coefficients now replaced by the modified ones Ea/s, taking care with the t 
and A dependences. This time the box at time vertex t corresponds to the sum 
—if^gpGapdt where now 

Ga/3 = (v-iVEV+(-i)'VEVEV + ...) E,p 

TTiVE)^/'^^) 

and this gives the required result. 

Step 5: Convergence of the Series 

What we have done so far has been to expand the Dyson series, determine 
the asymptotic limit of each diagram term (only the TC ones survived), to re- 
place the terms by their respective limits and to re-sum the series. To complete 
the proof, we need to establish that the series is absolutely and uniformly con- 
vergent. Fortunately we are able to extend proof for estimating these type of 
series exists for emission- absorption interactions ^01 to the general case. 

Let us start with the case where we have emission and absorption only in 
the interaction. The order must be even say, n — 2n2, as the vacuum diagrams 
consist of 71,2 pair contractions only. There will then be 2"2n2'- ^^'^^ diagrams 
with 2n2 vertices. A typical diagram, having n2 = 6 is sketched below: 



= -iEai3 - E^.V'i 



t2 h 



There exists a permutation a of the n — 2n2 time indices which re-orders to 
the diagram Dq (n) shown below 

The permutation is moreover unique if it has the induced ordering of the 
emission times. Not all permutations arise this way, the ones that do are termed 
admissible. We now consider an estimate of the n— th term in the Dyson series. 

Let E — max ||iJQ,^|| , then 



V •^'^"(*) ^ Admissible perms. '^'^W 

p "2 

E"" Y[\C{t2k-t2k-uX)\ 

J Bit) fc^i 
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where R (t) is the union of simplices {(in, • • • , ii) : t > tCT-i(n) > • • • > i(T-i(i) > O} 
over all admissible permutations Z. R{t) will be a subset of [0, i]^"^ and if we 
introduce variables t2k and S2fc — t2k — t2k-i for A: = 1, • • • , n2, it is easily seen 

max(i 1)"^ 

that the above is majorized by E'^^^C^^^ x , where C — maxV-' . 

This is the Pule inequality [30j and clearly gives the uniform absolute estimate 
required to sum the series. 

We now consider scattering, and constant, terms in the interaction. This 
time, the number of diagrams with n vertices will be given by the n-th Bell 
number Bn- To see why this is so, we recall that if we have a reservoir quanta 
created at a vertex, then perhaps scattered ,and finally reabsorbed then we can 
think of it as the same quantum and treat all the vertices it has been at as being 
linked. Each such diagram is then described by these subsets of linked vertices 
(we should also count the neutral vertices as these are singleton sets): in this 
way we have a one-to-one correspondence between the diagrams and partitions 
of vertices into non-empty subsets. The Bell numbers grow rapidly and have a 
complicated asymptotic behavior. The proliferation of diagrams is due mainly 
to the multiple scattering that now takes place. 

Let us consider a typical diagram. We shall assume that within the diagram 

there are ni singleton vertices [• • • , ■ ■ ■], n2 contraction pairs [• • • ....iT^"^.... 

• • • ], ns contraction triples [• • • ..../'^"^Y'^"^.... • • • ], etc. That is the diagram has 
a total of n = jnj vertices which are partitioned into m — connected 
subdiagrams. For instance, we might have an initial segment of a diagram 
looking like the following: 



There will exist a permutation a of the n vertices which will reorder the 
vertices so that we have the singletons first, then the pair contractions, then the 
triples, etc., so that we obtain a picture of the following type 

^3 triples pairs n\ singletons 

CYZv cx::\ ax c\ 

The permutation is again unique if we retain the induced ordering of the 
first emission times for each connected block. We now wish to find a uniform 
estimate for the n-th term in the Dyson series, we have 

V /" r[|C(A)| X "weights" 

2) "'A„(t) 

where the weights are the operator norms of various products of the type 
Ea^p^ ■ ■ ■ Each connected diagram of j > 2 vertices will typically have 
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one emission and one absorption, and j — 2 scattering vertices. The Pule argu- 
ment of rearranging the sum over diagrams into a single integral over a region 
R (t) of [0, t]" again applies and by similar reason we arrive at the upper bound, 
this time of the type 



E 



max(i,ir+"^+"=+- 



ni,n2,n3,--- ' ni!n2!n3! • • • 

Here the sum is restricted so that V • m, = n. An uniform estimate for the 
entire series is then given by removing this restriction; 

^ exp{j:^{Aj + B)nj} 
^{AB)- ^ ni!n2!n3!--- 

ni,n2,n3)--- 

where = ||VE|| and = Emax.{t, 1). Again we use the trick to convert a 
sum of products into a product of sums 

j n 

where we need < 1 to sum the geometric series - this however, is precisely 
our condition that ||VE|| < 1. 
Acknowledgment 

It is a great pleasure for the author to thank the staff at MabuchiLab for 

their kind hospitality during his visit there when part of this paper was writ- 
ten. Conversations with a Hideo Mabuchi, Gerard Milburn, P. Krishnaprasad, 
Matthew James, Howard Wiseman, Andrew Doherty, Ramon van Handel and 
Luc Bouten are gratefully acknowledged. 



References 

[1] C. Barnett, R. F. Streater, The Ito-Clifford Integrals, I. Wilde, J. Functl. 
Anal., 48, 172-212 (1982) 

[2] V. P. Belavkin, Theory of the control of observable quantum systems, Au- 
tomatica and Remote Control 44 (2) 178-188 (1983) 

[3] R. L. Hudson and K. R. Parthasarathy, Quantum ltd 's formula and stochas- 
tic evolutions. Commun. Math. Phys. 93, 301-323 (1984) 

[4] C. W. Gardiner and M. J. CoUett, Input and output in damped quantum 
systems: Quantum stochastic differential equations and the master equa- 
tion. Phys. Rev. A 31, 6, 3761 (1985) 

[5] E. B. Davies, Quantum Theory of Open Systems (Academic Press, London, 
1976) 



22 



[6] H. Spohn, Kinetic Equations from, Hamiltonian dynamics: Markovian lim- 
its, Rev. Mod. Phys. 52, 569 (1980) 

[7] H.M. Wiseman, G. J. Milburn, Quantum theory of optical feedback via 
homodyne detection, Phys. Rev. A 49, 4110 (1994) 

[8] R. van Handel, J. K. Stockton, H. Mabuchi, Feedback control of quantum 
state reduction, IEEE Trans. Automat. Control 50 (6), 768-780, (2005) 

[9] L. Bouten, S. Edwards, V. P. Belavkin, Bellman equations for optimal 
feedback control of qubit states, J. Phys. B, At. Mol. Opt. Phys. 38, 151- 

160(2005) 

[10] J. Gough, V.P. Belavkin, O.G. Smolyanov: Hamilton- J acobi- Bellman equa- 
tions for quantum filtering and control, J. Opt. B: Quantum Semiclass. Opt. 
7 S237-S244 (2005) 

[11] L. Accardi, A. Frigerio, Y. G. Lu, Weak coupling limit as a quantum func- 
tional central limit theorem, Commun. Math. Phys. 131, 537-570 (1990) 

[12] A. N. Chebotarev, Symmetric form of the Hudson-Parthasarathy equation. 
Mat. Zametki, 60, 5, 725-750 (1996) 

[13] J. Gough: Asymptotic stochastic transformations for non-linear quantum 
dynamical systems. Reports Math. Phys. 44, No. 3, 313-338 (1999) 

[14] V P Belavkin: Quantum stochastics, Dirac boundary value problem, and 
the ultra relativistic limit. Rep on Math. Phys. 46 (3): pp 359-381 (2000). 

[15] J. Gough, Quantum flows as Markovian limit of emission, absorption and 
scattering interactionsCommun. Math. Phys. 254, no. 2, 489-512 (2005) 

[16] M. Evans, R.L. Hudson Multidimensional Quantum Diffusions Quantum 
Probability III, Eds. L. Accardi and W. von Waldenfels, Lecture Notes in 
Mathematics 1303 69-88, (1988) 

[17] G. Lindblad: On the generators of completely positive semi-groups. Com- 
mun. Math. Phys 48, 119-130 (1976) 

[18] B. Oksendal, Stochastic Differential Equations, Universitext, Springer- 

Verlag, Berhn, (1992) 

[19] E. Wong, M. Zakai, On the relationship between ordinary and stochastic 
differential equations. Int.. J. Eng. Sci., 3, pp. 213-229 (1965) 

[20] C. W. Gardiner, Quantum Noise, Springer Series in Synergetics, (1991) 

[21] K. R. Parthasarathy, An Introduction to Quantum Stochastic Calculus, 
Monographs in Mathematics, (Birkhauser 1992) 

[22] E. N. Economou, Green's Functions in Quantum Physics, Springer Series 
in Solid-State Sciences, Volume 7 (1983) 



23 



[23] M. Reed & B. Simon, Methods of Modern Mathematical Physics, Vol 1., 
Theorem VIII. 3, (Academic press 1972) 

[24] J. Gough, Noncommumative Ito and Stratonovich noise and stochastic evo- 
lutions, Theoretical Math. Phys., 113, 276-284, (1997) 

[25] W. von Waldenfels, Symmetric differentiation and Hamiltonian of a quan- 
tum, stochastic process, Inf. Dim. Anal. & Quantum Prob. 8, No. 1, 73-116, 
(2005) 

[26] J. Gough, Quantum, White Noise and the Master Equation for Gaussian 
Reference States, Russ. Journ. Math. Phys., 10, No.2, 142-148, (2003) 

[27] W. H Louisell, Quantum Statistical Properties of Radiation, (Wiley Classics 
Library 1990) 

[28] R. Mattuck, A Guide to Feynman Diagrams in the Many-Body Problem, 
(Dover Science Books 1992) 

[29] J. Gough, A. Sobolev, Quantum Markovian approximations for Fermionic 
reservoirs. Inf. Dim. Anal. & Quantum Prob. 8, No. 3, 453-471, (2005) 

[30] J. V. Pule, The Bloch equations, Commun. Math. Phys. 38, 241-256 (1974) 



24 



